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Abstract. We review seven methods of measuring the absolute magnitude 
My of RR Lyrae stars in light of the Hipparcos mission and other recent 
developments. We focus on identifying possible systematic errors and rank 
the methods by relative immunity to such errors. For the three most robust 
methods, statistical parallax, trigonometric parallax, and cluster kinemat¬ 
ics, we find My (at [Fe/H]= -1.6) of 0.77 ± 0.13, 0.71 ± 0.15, 0.67 ± 0.10. 
These methods cluster consistently around 0.71 ±0.07. We find that Baade- 
Wesselink and theoretical models both yield a broad range of possible values 
(0.45-0.70 and 0.45-0.65) due to systematic uncertainties in the tempera¬ 
ture scale and input physics. Main-sequence fitting gives a much brighter 
My = 0.45 ± 0.04 but this may be due to a difference in the metallic- 
ity scales of the cluster giants and the calibrating subdwarfs. White-dwarf 
cooling-sequence fitting gives 0.67±0.13 and is potentially very robust, but 
at present is too new to be fully tested for systematics. If the three most 
robust methods are combined with Walker’s mean measurement for 6 LMC 
clusters, Vb,LMC = 18.98 ±0.03 at [Fe/H]=—1.9, then //lmc = 18.33 ±0.08. 


1. Introduction 

RR Lyrae stars are among the most popular local distance indicators. Smith 
(1995) describes in detail their general properties. One can measure the 
apparent magnitudes of RR Lyrae stars in a stellar system and infer their 
mean deredenned apparent magnitude Vq. If the mean absolute magnitude 
of RR Lyrae stars Mv{RR) at the system metallicity [Fe/H] is known, then 
the system’s distance modulus /r is: 


pi = Vo- Mv{RR). 


( 1 ) 
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In the following sections, we will assume that Vq can be measured accurately 
and concentrate on the My determinations. My can be calibrated through 
field or globular cluster stars. There are compelling arguments (e.g., Catelan 
1998) against two distinct, environment-dependent My scales (i.e., distance 
scales), and we will seek a universal absolute magnitude - metallicity My 
- [Fe/H] relation: 

My{RR) = a([Fe/H] + 1.6) + /?. (2) 

Most of the extragalactic distance scale is tied to the Large Magellanic 
Cloud (LMC). The My - [Fe/H] relation from (2) guarantees that Ulmc 
will be sensitive to the zero point (3 and insensitive to the slope a. Therefore, 
we do not discuss many a determinations (Ajhar et al. 1996; Fusi Pecci et al. 
1996; Kovacs &: Jurcsik 1996 etc.), but concentrate on methods yielding /3. 
Note that Vo,lmc ~ 19-0 (Walker 1992; Hazen & Nemec 1992; Reid h Freed¬ 
man 1994) and so, from equation (1) the division between the short and 
long distance scale occurs at /3 ~ 0.55, with fainter My corresponding to 
the short distance scale. Some methods determine My{RR) directly using 
their positions (trigonometric parallax), kinematics (statistical parallax) or 
pulsational properties (Baade-Wesselink). Others (cluster kinematics, main 
sequence and white dwarf fitting) establish globular cluster distances and 
then, in a second step, My{RR). We group methods according to their 
mathematical description. Our, rather incomplete, literature review serves 
mostly illustrative purposes, as we concentrate on the physical picture of 
the methods, emphasizing their strengths and weaknesses. 

2. Kinematic Methods 

The distance to an ensemble of stars can be measured by comparing their 
radial velocities (kms“^) and proper motions (s“^) under the assumption 
that these are due to (statistically) identical physical velocities. In effect, 
one fixes the distance so that either the first moment of the population (bulk 
motion) as determined from the proper motions is equal to first moment 
as determined from the radial velocities, or so that the second moments 
(dispersions) are equal, or both. 

The great beauty of these methods is that the basic measurements are of 
dimensionless quantities (redshift for radial velocities and positions on the 
sky for proper motions) and therefore no assumptions about the distance 
scale enter the determinations. The major systematic uncertainty (aside 
from concerns about the quality of the data) is that it may be difficult to 
verify that the radial velocities and proper motions in fact arise from the 
same physical velocity distribution. 

One may show that if the measurement errors are small compared to 
the intrinsic dispersion of the population, then the fractional distance error 
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from equating the first moments is 5r]/rf = l/[(yliAI)^/^Ksin0] where N 
is the total number of stars, 0 is the angular size of the system, Ai is 
a geometrical factor of order unity, and k is the ratio of bulk motion to 
velocity dispersion. Similarly the error from equating the second moments 
is drj/r] = l/(A 2 iV)^/^ where A 2 is another factor of order unity. Hence, by 
combining the two methods, the error is 


N[Ai {k sin Oy + A 2 ] 

For nearby open clusters like the Hyades, k r\j 10^ and 0 ~ 10 so 
the first term in the denominator dominates, and the distance is deter¬ 
mined from the first moments. This is sometimes called the “moving clus¬ 
ter” method. However, there are no RR Lyrae stars in open clusters. For 
globular clusters, k ~ 10 ^ and 6 ~ 10 “^, so the second term dominates. 
Hence, kinematic distances to globular clusters are determined by match¬ 
ing velocity dispersions. For field RR Lyrae stars in the solar neighborhood, 
K ~ 2.1 and sin0 ~ 1, so both the first and second moments are used to 
determine the RR Lyrae distance scale. This is called “statistical paral¬ 
lax”. Statistical parallax and kinematic distance to globular clusters have 
different sources of systematic error, so we treat them separately. 

2.1. STATISTICAL PARALLAX 

Statistical parallax is reviewed thoroughly by Layden elsewhere in this vol¬ 
ume, so we give only a brief overview here. In contrast to the case of 
globular clusters, local field RR Lyrae stars are not at a common dis¬ 
tance. Hence, before their radial velocities and proper motions can be 
compared, the stars must be put on a common scale by taking account 
of their dereddened apparent magnitudes. That is, one can measure the 
9 parameters describing the velocity ellipsoid (3 components of bulk mo¬ 
tion plus 6 independent components of the velocity-dispersion tensor) from 
radial velocities alone. On the other hand, if one assumes some arbitrary 
Mv{RR), for the RR Lyrae stars, then one can infer their distances from 
their measured apparent magnitudes and estimated extinctions. The dis¬ 
tances and proper motions yield the transverse velocities, and from these 
one can again estimate the 9 parameters of the velocity ellipsoid. One could 
then adjust the assumed My so that the velocity ellipsoid from proper 
motions matched the velocity ellipsoid from radial velocities as closely as 
possible. In practice, one fits for all 10 parameters (9 for the velocities 
plus My) simultaneously using maximum likelihood. The maximum like¬ 
lihood approach was pioneered by Clube & Dawe (1980) and was subse¬ 
quently applied by Hawley et al. (1986) to the then best-available data 
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set of 142 RRab stars. They obtained My = 0.76 ± 0.14. Note that the 
error is very close to the theoretical minimum (for no measurement er¬ 
rors) aMv = (5/ln 10)[(2A^/9)(6-|-= 0.12 mag (Popowski & Gould 
1998a). This is because the measurement errors do not contribute signifi¬ 
cantly if they are substantially below the velocity dispersion. 

As we discuss below, this result is fainter than virtually all other es¬ 
timates and much fainter than some. Since the method itself appears ex¬ 
tremely robust, several workers have invested substantial effort to obtain 
the most reliable input data and to investigate whether any unrecognized ef¬ 
fects could be leading to systematic errors. Layden (1994) and Layden et al. 
(1996) put all existing data on a homogeneous system. They incorporated 
the proper motions from the new Lick NPMl (Klemola, Hanson, & Jones 
1993) survey which has smaller (and just as important, better-understood) 
errors than previous proper-motion studies. They used the Lick catalog to 
calibrate the errors of the non-Lick proper motions. They found that the 
diverse photometry sources were offset from one another and put these on a 
common zero point. Likewise, they put all the extinctions on the Burstein 
& Heiles (1984) system. They found My = 0.71 ± 0.12 for 162 “halo-3” 
RRab stars with mean metallicity ([Fe/H]) = —1.61. Popowski & Gould 
(1998a) developed a formulation of maximum likelihood which permitted 
both new analytic investigations and much more vigorous Monte Carlo in¬ 
vestigations of possible systematic effects. These included a possible change 
in the velocity ellipsoid with distance from the Galactic plane, rotation of 
the local coordinate frames relative to the Sun’s frame, and effects due to 
the severe non-Gaussianity of the velocity distribution. However, these ef¬ 
fects all proved negligible. The most important previously overlooked effect 
that they found was Malmquist bias, and they obtained My = 0.76 ± 0.12 
for the 162-star Layden et al. (1996) sample. 

Hipparcos has had two major impacts on the RR Lyrae statistical par¬ 
allax determination. One is, of course, new and more precise proper-motion 
measurements. The other, more surprisingly and indeed more importantly, 
is better photometry. Fernley et al. (1998a) fit Hipparcos light curves to 
obtain new photometry for most RR Lyrae stars in the Hipparcos cata¬ 
log. They combined these with Hipparcos proper motions and obtained 
My = 0.77±0.17 for 84 “halo” ([Fe/H]< —1.3) stars with ([Fe/H]) = —1.66, 
including 69 RRab’s and 15 RRc’s. Fernley et al. (1998a) also used high- 
precision ground-based photometry to show that their Hipparcos-based 
mean magnitudes were correct with very small (< 0.02 mag) scatter. (Tsu- 
jimoto, Miyamoto, & Yoshii 1998 conducted a similar study of 99 Hipparcos 
“halo” stars and found My = 0.69 ± 0.10 at ([Fe/H]) = —1.58. However, 
since their quoted error is a factor 0.7 below the theoretical minimum, we 
conclude that their analysis is incorrect.) Popowski & Gould (1998b) used 
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the Hipparcos proper motions to check earlier catalogs and found that only 
the Lick catalog is of sufficiently high quality to use. Gould & Popowski 
(1998) used Fernley et ah’s (1998a) Hipparcos-based mean magnitudes to 
check Layden et al.’s (1996) systematization of previous heterogeneous pho¬ 
tometry, and found that it was 0.06 mag too bright. They also incorporated 
the new extinction map of Schlegel, Finkbeiner & Davis (1998) based on 
COBE/IRAS measurements of dust emission, and they eliminated a num¬ 
ber of stars with questionable extinctions and proper motions to obtain 
My = 0.77 ± 0.13 for 147 “halo-3” RRab stars with ([Fe/H]) = —1.60. 

Layden et al. (1996) noted that the velocity ellipsoid of their solution is 
in good agreement with that of Beers & Sommer-Larsen (1995) for metal- 
poor stars and took this as independent confirmation of the correctness of 
their results. Popowski & Gould (1998b) and Gould &: Popowski (1998) 
directly incorporated the Beers & Sommer-Larsen stars into the analysis 
and obtained My = 0.80±0.11 at ([Fe/H]) = —1.71 for a combined sample 
based on 149 RRab Lyrae stars and 716 non-RR Lyrae stars. 

At this point, essentially all systematic errors have been removed from 
the RR Lyrae statistical parallax determination. The statistical errors in 
the above two solutions (My = 0.77 ± 0.13 at ([Fe/H]) = —1.60 or My = 
0.80 ± 0.11 at ([Fe/H]) = —1.71) should therefore be taken at face value. 

2.2. KINEMATIC CLUSTER DISTANCES 

If the distance to a cluster is known, the My{RR) can be determined by sub¬ 
tracting the distance modulus from the dereddened apparent magnitude of 
RR Lyrae stars in the cluster, or more generally from the height of the zero 
age horizontal branch (ZAHB) at the color of the instability strip. Cluster 
distances can be determined kinematically by comparing the dispersions of 
the radial velocities and proper motions. The principle is similar to statisti¬ 
cal parallax, but there are two major practical differences. First, the proper 
motions are much smaller for clusters (~ 10kms“^/10kpc = 0.2masyr“^) 
than for nearby field RR Lyrae stars (~ 200kms“^/2kpc = 20masyr“^). 
Second, clusters are seen in projection, so the 3-space position of individ¬ 
ual stars is unknown. This introduces additional systematic effects that are 
difficult to fully take into account. 

Cudworth (1979) made the first such measurement, finding a distance 
d = 9.6 ± 2.6 kpc for M3 by comparing the proper-motion dispersion of 
71 stars, = 0.094 ± 0.021 mas yr“^ with the radial-velocity dispersion 
from Gunn & Griffin (1979), ar = 4.3 ± 0.7kms“^. The measured rms dis¬ 
persion was am = 0.183 mas yr“^ while the “internal” rms measurement 
error was e = 0.130masyr“^. Cudworth took the “external” error to be 
/e and estimated / = 1.20 ± 0.12. Note that the measurement errors were 
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actually larger than the cluster dispersion = ci^ — (/e)^ and that the 
correction factor / is therefore very important. This factor had earlier been 
measured as / = 1.1 ± 0.1 by Cudworth & Monet (1979) by comparing 
reductions of two different plate sets of M13, one weak and one deep. How¬ 
ever, Cudworth (1979) argued that this was only a lower limit on / since 
the two sets were taken with the same Yerkes telescope, over the same time 
interval, ca. 1900 to ca. 1975. We note that a perfectly plausible / ^ 1 
would yield ci^ = 0.129 ± 0.023mas yr“^ and hence a cluster distance of 
d = 7.1 ± 1.7kpc. 

Cudworth (1979) assumed that M3 is isotropic because it looks circu¬ 
larly symmetric on the sky and rotation is barely detectable if at all in the 
radial velocities. Hence, he simply divided the radial-velocity dispersion by 
the proper-motion to obtain the distance. However, if the cluster rotation 
axis were along the line of sight, these observational characteristics would 
remain in tact, but the radial dispersion would no longer be representa¬ 
tive of the tangential dispersion. (Transverse rotation cannot be measured 
from the proper-motion data because there are too few foreground stars 
and their dispersion is too large to form a stable framework.) Thus, M3 
exemplifies both the promise and the problems of this technique. 

Gunn & Griffin (1979) introduced a major advance by incorporating 
Mitchie-King models (Mitchie 1963) into the analysis, and Lupton, Gunn, 
& Griffin (1987) built on their experience when they measured the distance 
to M13. The Mitchie-King models are constrained by the light profile, the 
radial-velocity map, and the mass function. Such models can take account of 
rotation although the degeneracy noted in the previous paragraph remains. 
M13 is closer than M3 (6.5 vs. 9.6 kpc), the velocity dispersion is larger 
(6.5 vs. 4.3 kms“^), and the number of stars is larger (268 vs. 71). Hence 
the distance errors are much smaller (8% vs. 27%). Lupton et al. (1987) do 
not say whether they incorporated Gudworth & Monet’s (1979) estimate 
of / = 1.1, or whether they used / = 1. In this case, the difference in the 
final result is 4%, or 0.08 mag in the height of the horizontal branch (HB). 
Lupton et al. (1987) determined the distance only for their best-fit model 
of the cluster and did not obtain distances for a range of acceptable models 
which would have allowed them to evaluate the systematic uncertainty due 
to their modeling procedure. 

Peterson & Gudworth (1994) measured the distance to M22 by equating 
the minor axis proper-motion dispersion with the radial-velocity dispersion 
and found d = 2.6 ± 0.3 kpc. That is, they assumed that the axis of the 
observed rotation is in the plane of the sky and, rather than attempt to 
model the effect of this rotation on the observed proper motions in the 
symmetry plane, simply ignored half the proper-motion data. They do not 
say whether they incorporated an / factor, but for this case the difference 
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between / = 1 and / = 1.1 is only 0.03 mag. 

Rees (1996) has applied Gunn & Griffin’s (1979) Mitchie-King tech¬ 
nique to 8 clusters (47 Tuc, M5, M4, M92, N6397, M22, M15, and M2) 
and combined the results with the distances obtained for M3 and M13 
by Cudworth (1979) and Lupton, Gunn &: Griffin (1987). Subsequently, 
Rees (1998 priv. comm.) rereduced the M22 data and derived a distance, 
d = 11.26 ± 1.31 kpc. He also came to the conclusion that the N6397 proper 
motions were not reliable enough to use. Two of the remaining 9 clusters 
have features which are not well represented by Mitchie-King models: M15 
has a central cusp, and 47 Tuc has differential rotation. If these two clusters 
are eliminated, the remaining seven have a mean zero-age HB (ZAHB) of 
My = 0.62 ± 0.10 mag at ([Fe/H]) = —1.6. (Rees quotes an uncertainty of 
0.05 mag, evidently having renormalized the errors because = 1-66 for 6 
degrees of freedom. However, since the errors are basically due to counting 
statistics - and so cannot have been overestimated - the errors should not 
be renormalized.) If M15 and 47 Tuc are re-included. My = 0.58 ± 0.08. 

The details of Rees’s (1996) work are not yet publicly available, so it is 
impossible to give a complete assessment of the systematic errors. Glearly, 
if the proper motion errors have been systematically underestimated (as 
suspected by Gudworth & Monet 1979 and Cudworth 1979), then the lu¬ 
minosity of the ZAHB has been underestimated for each of the clusters. 
This is a much smaller effect for nearer clusters which have smaller over¬ 
all errors and so more statistical weight. Alternatively, if there is no error 
underestimate but an underestimate has been corrected for, then the lumi¬ 
nosity of the ZAHB has been overestimated. To date, no one has published 
a set of distances for a set of acceptable Mitchie-King models, so there is no 
way to assess the systematic error due to uncertainty in the cluster geom¬ 
etry. At least until such tests are conducted, clusters that are inconsistent 
with Mitchie-King models (i.e. M15 and 47 Tuc) should be excluded. 

Finally, there is bias introduced by binaries which increase the radial- 
velocity dispersion but not the proper-motion dispersion (measured on 100 
yr time scales). The correction for this previously unrecognized bias makes 
My fainter by 5My = (2.5/In 10)/((7V)5rb)^)/3u'^ niag, where / is the bi¬ 
nary fraction, a is the 1 -dimensional dispersion of the cluster, V(,rb is the 
mean orbital velocity of the binaries, 7 = Ms/{Ms + Mp), and Mp and Mg 
are the primary and secondary masses. Hut et al. (1992) conclude (primar¬ 
ily based on Pryor et al. 1989) that for periods of 0.2 to 20 years, the binary 
fraction of globular clusters is consistent with that of G stars as measured 
by Duquennoy & Mayor (1991), i.e., 8 %. Applying the above formula to 
the Duquennoy & Mayor (1991) distribution, we find 6 My = 0.05. Hence, 
our best estimate for this method is My = 0.67 ± 0.10 at ([Fe/H]) = —1.6. 
We caution that this error is only statistical and that the systematics are 
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3. Trigonometric Parallax 

Prior to Hipparcos, it was not possible to make a useful estimate of the 
Mv{RR) using trigonometric parallax. Hipparcos measured only one RR 
Lyrae star with reasonable precision, RR Lyrae itself. This yields My = 
0.78±0.29 at [Fe/H]= —1.39 (Fernley et ah 1998a), still not precise enough 
to discriminate among various other estimates. 

However, Gratton (1998) has extended this approach by considering all 
HB stars (not just RR Lyrae) in a magnitude limited sample, Vq < 9. This 
criterion yields 22 stars, including 10 blue HB stars, 3 RR Lyrae stars, and 9 
red HB stars. Gratton (1998) obtains My = 0.69 ±0.10 at (Fe/H) = —1.41 
if all 22 stars are included, or My = 0.60±0.12 at (Fe/H) = —1.51 if one red 
HB star that he suspects of being a giant (not HB) star is eliminated. We 
first review and correct Gratton’s statistical procedure, and then discuss 
more carefully the problem of selection. 

Gratton (1998) begins by determining the shape of the HB from observa¬ 
tions of M5. This yields 6Mv{B — V), the difference between the magnitude 
at the instability strip and the magnitude of a star at a particular B — V 
color. His model of the HB is therefore characterized by a single parameter. 
My at the instability strip. From the observed magnitude and color and 
the inferred reddening, he is then able to predict the parallax as a function 
of this parameter, 

7r*(My) = lo0.2{[Mv-5Mv{B-V)]-[V-3.1E{B-V)]+10}^g^^ ^ 4 ^ 


He then forms the average of the difference between this quantity and the 
observed parallax, tt, weighted by the observational errors, and determines 
My by setting this quantity to zero 




( 5 ) 


Now, the proper procedure would be to form x^iMy) = i^v) — 

7 ri)‘^la‘^^ and then to minimize it with respect My. Formally, 


^ dx^jMy) 

~ dMy 


In 10 TT* (My) — TTj 


2.5 


i 


(J, 


<{My). 


7r,2 


( 6 ) 


Comparing equations (5) and (6), we see that Gratton (1998) has in effect 
weighted the terms by their distance, which of course gives higher weight to 
the most poorly determined values. Employing equation (6) rather than (5), 
we find My = 0.75 ± 0.09 and My = 0.64 ± 0.12 for the two cases. 
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The red HB star that Gratton (1998) eliminated is suspicious only be¬ 
cause it is significantly fainter than the HB. Of course, since its paral¬ 
lax is only 1.3 a larger than predicted, this could be a normal statisti¬ 
cal fluctuation. There does not appear to be any way to guard against 
giant-star contamination of the red HB sample without biasing the sample 
by eliminating stars with above average, but normally distributed, paral¬ 
laxes. Therefore the most prudent procedure is to eliminate all red HB 
stars. We then find My = 0.68 ± 0.14 at (Fe/H) = —1.62. Finally, we 
note that account must be taken of the intrinsic scatter of the population. 
For homogeneous populations (e.g. globular clusters) the scatter in My is 
aMv ~ 0-08 mag. The upper limit for the scatter in a general field popula¬ 
tion is aMv ~ 0.17 mag, but a more plausible estimate is aMy ~ mag 
(Popowski & Gould 1998a). We adopt that here and find My = 0.68 ±0.15. 
This scatter introduces a Malmquist bias. One may show that if the par¬ 
allax errors scale with flux as a-,^ oc F~'^ , then the Malmquist bias is 
5My = (lnl0/5)(l - i^)alj^ 0.01 mag, where u ~ 0.2 (appropriate for 
6.5<F<9 for Hipparcos). Our final estimate is therefore My = 0.69 ± 0.15 
at (Fe/H) = —1.62, where the metallicity is on the Gratton scale, corre¬ 
sponding to My = 0.72±0.15 at (Fe/H) = —1.6 on the Zinn & West (1984) 
scale. 


4. Baade-Wesselink Method 

The goal of the Baade-Wesselink method is to determine the equilibrium ra¬ 
dius of a pulsating star through the combined analysis of its radial-velocity 
and light curves. The equilibrium radius and effective temperature, , 
inferred from multicolor photometry allows one to determine the luminos¬ 
ity, L. The radial-velocity curve reflects the absolute change in radius R 
of a star, 5R, due to its contraction or expansion with velocity Vrad{t)- 
Gonsequently, during the time interval 6t: 


5R{t, t + 5t) = Vrad{t, t + St) * St, 


( 7 ) 


where Vrad{t-, t + St) is the mean radial velocity in a time interval (t, t + 5t). 
The flux of a spherical blackbody at distance d is given by; 

F - - 55 -■ ( 8 ) 

where as is a constant. Consequently, one can determine R by: 


R = 25R 


( T^siR)F{R + SR) 
[t^^{R + SR)F{R) 



-1 


(9) 
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In equation (9), SR is known from radial velocity curve [eq. (7)], the fluxes 
are observables, and the temperature at each phase is estimated from colors 
based on stellar atmospheric models (e.g., Kurucz 1992). 

Most of the results obtained before 1994 for field RR Lyrae stars (Liu & 
Janes 1990; Jones et al. 1992; Cacciari, Clementini & Fernley 1992; Skillen 
et al. 1993) can be summarized by one equation; 

My = 0.70+ 0.21 ([Fe/H]+ 1.6), (10) 


resulting in ;Ulmc ~ 18.3. 

The Baade-Wesselink method is more involved and model-dependent 
than geometric and kinematic determinations. The simple physical picture 
of the method presented above gives an immediate insight into its possible 
systematics. First, Vrad in equation (7) is not the mean radial velocity 
that is measured from spectral lines. Vrad would be equal to the measured 
velocity if the stellar photosphere were moving straight toward or away 
from the observer. In fact it is moving radially inward or outward from the 
center of a star and its direction forms some angle, 9, with the observer- 
star line of sight, which reduces the observed radial velocity by a factor 
of cos(0). Additionally, a contribution of a given stellar surface patch to 
the line intensity depends on limb darkening of the star and the relative 
depths at which a given line and the stellar continuum form. All these effects 
are usually parametrized by a factor “p” which is supposed to convert the 
measured radial velocity to a true one (e.g., Fernley 1994). Getting (1935) 
obtains p = 1.41 analytically using the Milne-Eddington limb darkening 
model. More sophisticated investigations using model atmospheres usually 
produce slightly lower values in the range 1.28 - 1.39. The factor p is a 
function of pulsation phase, but Fernley et al. (1989) showed for X Ari that 
using a varying value of p changes the derived R by less than 1%. The factor 
p also increases sharply with decreasing line strength (Karp 1975). The 
uncertainty in p does not influence the slope of the My - [Fe/H] relation 
but does have a substantial impact on the zero point because F oc p^ 
[eqs. (9), (7) into (8)]. For small changes in p, the zero-point shift is 


A(3 = 


5 Ap 

In 10 p 


( 11 ) 


Fernley (1994) argues that the original values of RR Lyrae p coefficients 
with a mean of ~ 1.32 (Skillen et al. 1993) should all be replaced with a 
universal value of p = 1.38. Equation (11) then predicts A/3 = 0.09, but 
more detailed analysis (Eernley 1994) gives A/3 = 0.07, leading to 


My = 0.63 + 0.21 ([Fe/H] + 1.6) . 


(12) 
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Second, the unambiguous determination of the effective temperature 
from stellar colors is still problematic. Temperature differences are likely 
to affect both the slope and zero point of the My - [Fe/H] relation. Tgff as 
obtained from optical colors (A £800 nm) are higher by 200-300 K than 

those inferred from (V — K) colors. Most observers prefer (V — K) and 
usually argue that the existence of {U — B) excess in some pulsational 
phases is an indication of atmospheric shocks. Therefore, optical colors are 
likely to give an incorrect Tgfj . However, McNamara (1997) claims that 
Preston & Paczynski (1964) find no correlation between ultraviolet excess 
and strength of hydrogen emission lines suggesting that some other phe¬ 
nomenon is responsible for this effect. He also argues that surface gravities 
as inferred from optical colors match much better the expectations from 
mass and radius values for RR Lyrae stars, whereas the surface gravities 
derived from (V — K) are too small. McNamara (1997) derives an My - 
[Fe/H] relation with a rather steep slope and bright zero-point; 


My = 0.50-^0.29 ([Fe/H]-M.6). (13) 


Which set of colors should be used to obtain the My - [Fe/H] relation re¬ 
mains an open question. Feast (1997) argues that even “traditional” Baade- 
Wesselink analyses produce steep slopes. He claimes that My should be 
treated as the independent variable rather than [Fe/H] . This is incorrect. 
If a regression is done, the quantity with smaller relative errors should be the 
independent variable. Maximum likelihood, which treats both errors equally 
is the most rigorous approach and reproduces the shallow slope (Fernley et 
ah 1998b). In addition. Feast’s (1997) claim that we do not see the brightest 
RR Lyrae stars due to a bias of the same sort of Malmquist bias is unlikely. 
A more natural explanation for the lack of stars with My < 0.5 is that they 
do not exist; the My - [Fe/H] relation may simply flatten at low [Fe/H] as 
predicted by some theoretical models (e.g. Caputo 1997). 

Third, RR Lyrae stars are not perfect blackbodies and bolometric cor¬ 
rections must be applied. Even though RR Lyrae stars are rather hot and 
metal poor, so departures from blackbody should be modest, the adopted 
bolometric corrections actually differ by up to 0.1 mag (Sandage 1993 vs. 
Kurucz 1992). 

Similar investigations have been performed for globular cluster RR Lyrae 
stars. Representative results obtained by Storm, Carney & Latham (1994), 
namely My = 0.66 at [Fe/H] = —1.4 for M5 and My = 0.45 at [Fe/H] = 
—2.24 for M92, are slightly brighter than values predicted by equation (10). 



12 


PIOTR POPOWSKI AND ANDREW GOULD 


5. Theoretical Models of Horizontal Branch 

One can study the My — [Fe/H] relation using a purely theoretical ap¬ 
proach by constructing mock HBs in the Tefj — L plane and evolving them 
in time. The results depend sensitively on the input physics (especially the 
opacities and equation of state: Castellani, Brocato & Persimoni 1996), He 
abundance, and assumed core mass. The treatment of convection is also one 
of the most essential ingredients. To compare theoretical HBs to HBs ob¬ 
served in globular clusters one must transform from the theoretical Tefj — L 
plane to the My - color plane. This step requires accurate atmosphere mod¬ 
els (e.g. Kurucz 1992). Then one can investigate HB morphology, how HB 
appearance depends on different parameters (e.g., age) etc. One can also 
establish how Mv'(ZAHB) at the instability strip scales with the [Fe/H] of 
a given stellar ensemble and how this relation depends on the evolution¬ 
ary stage of the population. With some (even conservative) observational 
constraints taken into account, the model predictions become quite robust. 
Another approach involves the analysis of the pulsational properties of RR 
Lyrae stars in the spirit of the pioneering work by van Albada &: Baker 
(1971). Pulsational models predict the relations between Teff, [Fe/H], lumi¬ 
nosity, L, period, P, and mass, M, of the stars; e.g., 

log P = 11.497 -F 0.84 log L - 0.68 log M - 3.48 log (14) 

They also enable one to investigate the location of blue and red, fundamen¬ 
tal and first overtone instability strip boundaries (Bono et al. 1995). Which 
regions of the instability strip are populated by RRab and RRc stars corre¬ 
late with the Oosterhoff type of the cluster and gives insight into the source 
of the period-shift effect (Sandage 1981,1982; Bono et al. 1995) which has 
been used to determine the My — [Fe/H] relation. The most robust con¬ 
straints on the My - [Fe/H] relation come from combining the pulsational 
and evolutionary scenarios. 

The predictions of theoretical evolutionary models are sensitive to every 
change in input physics. For example, recent improvements in opacities 
and the equation of state led to an increase in HB luminosities by ~ 0.1 
mag. Gratton et al. (1997) report an unpublished My - [Fe/H] relation by 
VandenBerg et al. (1998): 

My{ZAHB) = 0.19([Fe/H] + 1.6) + 0.66 (15) 

Salaris et al. (1997) obtain a slightly brighter result with a similar slope: 

My{ZAHB) = 0.21([Fe/H] + 1.6) + 0.57 (16) 

However HB luminosities are very sensitive to the He core mass. Mazzittelli 
et al. (1995) suggested that core masses should be increased by O.OIMq 
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relative to older estimates. Models with heavier cores have slightly steeper 
and brighter HBs, e.g. Caloi et al.’s (1997) non-linear My - [Fe/H] relation 
which can be approximately characterized by 

MviZAHB) = 0.26([Fe/H] + 1.6) + 0.49 (17) 

Caputo et al.’s (1993) and slightly adjusted Caputo (1997) HB models 
predict a break in slope which occurs near [Fe/H] = —1.3 if RR Lyrae stars 
have scaled solar abundance of a elements, and near [Fe/H] = —1.6 if a 
elements are enriched by a factor of 2-3. These models predict a steep slope 
of 0.3 for metal-rich RR Lyrae stars and a shallow slope of 0.19 for metal- 
poor stars. The zero point favors the short distance scale with ulmc ~ 18.4. 
Caputo (1997) checks her synthetic HB evolutionary models proving that 
they satisfy the constraints produced by convective pulsating models with 
a range of masses and He abundance Y = 0.24. Such models allow one 
to predict the borders of the RR Lyrae distribution in the P - My plane. 
Moreover, Caputo (1997) concludes that adopting the “long distance scale” 
with a steep slope for all metallicities (My = 0.3[Fe/H] -|- 0.94) results in 
RRc stars falling in the hot stable region, i.e. outside the instability strip. 
On the other hand, the calibrations with a shallow slope are in agreement 
with observations only if the zero point suggested by Fernley (1994) is 
adopted instead of a fainter value. 

Lee, Demarque &: Zinn (1990) and Caputo et al. (1993) analyzed HB 
evolutionary tracks showing explicitly how RR Lyrae luminosity depends 
on evolutionary stage and that blue, less-massive HB stars cross instability 
strip at significantly higher L than the ZAHB. Therefore the average RR 
Lyrae luminosity will depend on which part of the HB they come from 
and so what is the HB morphology. Particular My — [Fe/H] relations will 
therefore in general depend on the HB morphology of the cluster. If the 
cluster contains B blue HB stars, V RR Lyrae stars, and R red HB stars, 
then its morphology index (Lee 1989), (B — R)/(B-l-V + R), is a number in 
the interval (—1,1). Morphologies below -0.8 (47 Tuc) and above 0.8 (M92) 
indicate very red and very blue HBs and, according to Caputo (1997), 
are likely to complicate the picture. However, even if one restricts oneself 
to intermediate morphologies, the uncertainties of HB models (at least 0.1 
mag/dex in slope and 0.2 mag in zero point) prohibit any precise conclusions 
about the theoretically-preferred My - [Fe/H] relation. 

There is also another theoretical path leading to a determination of My. 
It progresses through the analysis of double mode RR Lyrae stars which 
is rather complex and is beyond the scope of this paper. Representative 
results (Alcock et al. 1997, Kovacs &: Walker 1998) tend to give a bright 
My. 
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6. Main or White-Dwarf Sequence Fitting to Globular Clusters 

The main sequence (MS) or white-dwarf (WD) cooling-sequence fitting 
techniques give the distance to a cluster by matching the cluster sequence 
of subdwarfs or WDs to a corresponding sequence in My - dereddened 
color plane as determined from calibrating stars in the solar neighborhood. 
The My of the HB can be determined from the equation: 


Mv{RR) 


V{RR) - Av{RR) - 
V{RR) - 


V{CS) - Mv{LC) - Av{CS) + ^ 


dMy 

d'ji 


V{CS) - Mv{LC) + V ^(57, 

^ d'ji 




( 18 ) 


where V{RR) and V{CS) are the apparent magnitudes of RR Lyrae stars 
in a cluster and the cluster sequence, respectively, Ay{RR) and Ay{CS) 
are the visual extinctions toward them, and My{LC) is the My of the lo¬ 
cal calibrators at a given dereddened color. Finally, J2i corrects for 

the fact that for cluster stars it is not always possible to find local cali¬ 
brators with the same characteristics (parameter 7 *) or just for parameter 
uncertainty of order of d'ji. For MS fitting the most essential parameters 
are metallicity and color (e.g., {B — V)), for WD fitting — mass and color. 
Note that if the cluster sequence stars experience the same extinction as 
RR Lyrae star, then in equation (18), Ay{RR) and Ay{CS) cancel out. In 
practice it is enough if MS or WD sequence stars are drawn from the same 
region of the cluster as RR Lyrae stars. There are a few crucial criteria that 
must be met to make the fitting method work properly: 

1. accurate distances to the local calibrators must be known. Only when 
one knows the distances is it possible to construct MS or WD curves 
in the My - dereddened color plane. In principle one should also know 
the reddenings to individual stars, but for the calibrating stars, most 
of which sit in the Local Bubble, this is not a major concern. 

2 . if the results of the sequence fitting are sensitive to some parameter, 
numerous stars in different ranges of this parameter should be avail¬ 
able. Metallicity is a key parameter for the MS fitting because the My 
of the MS at the given color is a function of metallicity, especially at 
the high metallicity end. One wants to avoid using theoretical models 
to shift MS defined at one metallicity to get the My for another. 

3. one must measure the parameters of the cluster on the same scale as 
applied to the local calibrators, e.g., the metallicity of a cluster should 
be measured on the same scale as the one used for local subdwarfs. It is 
necessary to match the appropriate local sequence to the one observed 
in a cluster. One possible cause of a systematic error in the MS fitting 
distances is that the metallicities of the local subdwarfs might be on 
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a different scale from those of the clusters (determined from giants). 
Even if the procedure to obtain each type of metallicity is completely 
uniform, the conclusions depend on the atmospheric models, which are 
completely different for subdwarfs and giants. Specifically, if the subd¬ 
warf metallicities were too low (or the giant metallicities too high) then 
the cluster distance and the HB luminosities would be overestimated. 
King et al. (1998) found intriguing evidence of a possible misalign¬ 
ment of this sort. They measured the metallicities of M92 subgiants 
(unfortunately not subdwarfs, but with higher gravity than giants), 
and obtained metallicities up to half a dex lower than those of M92 
giants. Pinnsoneault (1998, priv. comm.) finds a metallicity of M5 from 
BVI colors of MS stars. His value is ~ 0.3 dex more metal-poor than 
Gratton et al.’s (1997) giant-based determination. 

4. reddening corrections to the color must be known very accurately. Red¬ 
dening changes the conclusions about the My of MS or WD stars in 
two major ways. First, through the extinction Ay, it directly dims the 
light of the observed stars shifting the cluster sequence vertically on 
the color magnitude diagram. However, this effect almost entirely can¬ 
cels the Ay correction when one determines the My of the HB — both 
types of stars are extincted by approximately the same amount. And 
second, uncertainty in color, 5{X — Y), results in the shift in My of 
order of SMy ~ d{^-Y) '5{X — Y). This effect is particularly dangerous 
if the My - color relation is very steep (e.g., for [B — V) color of MS 
stars, 5My ~ —5 • 5{B — H)). 

6.1. MAIN SEQUENCE FITTING 

We stress at the outset that the MS fitting technique is very sound when 
applied to nearby open clusters. For a given metallicity there are plenty 
of local calibrating dwarfs that define well constrained sequences. Because 
they are numerous there is no problem with choosing a sample with very 
accurate parallaxes and so small errors in My. On the cluster side - be¬ 
cause open clusters are close, the spectroscopy of the MS stars is straight 
forward. Therefore, the clusters and local calibrators are easily put on the 
same metallicity scale. Globular clusters pose much more severe problems. 
Their metallicities are generally very low and, consequently, there are few 
local subdwarfs that can define metal-weak MS’s. The local subsample is 
collected from a larger region, which translates into larger parallax er¬ 
rors. Therefore some of the distances to individual stars require a sub¬ 
stantial Lutz-Kelker type correction (Lutz & Kelker 1973; Hanson 1979; 
Smith 1987). Additionally, globular cluster are far away and their MS’s 
are too faint for high-resolution spectroscopy. Therefore the cluster metal- 
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licities come from giants and not from subdwarfs. Reddening to globular 
clusters is usually measured from integrated photometry (e.g., Zinn 1980; 
Reed, Hesser, &: Shawl 1988), and potentially more accurate Stromgren 
photometry is generally available only for southern clusters. Moreover, for 
unknown reasons, Stromgren-photometry-based reddenings are systemati¬ 
cally smaller than other determinations. As we mentioned above, reddening 
uncertainties can dramatically affect the distance modulus. 

Reid (1997) uses 15 subdwarfs with Hipparcos parallaxes accurate to at 
least 12% to dehne local subdwarf sequences of known My. He separates the 
calibrating stars into two metallicity bins; intermediate with [Fe/H] ~ —1.4 
on the Zinn & West (1984) scale and metal-poor with [Fe/H] ~ —2.1. That 
is, Reid (1997) does not apply any theoretical transformations to obtain 
the My of MS stars at a given color. He reports that observational se¬ 
quences are much less sensitive to metallicity than predicted by theoretical 
isochrones. However, this effect may be spurious because Reid (1997) uses 
Carney (1994) metallicities that are based on low-resolution spectra and 
therefore are likely to show more metallicity spread than is really present. 
He derives new distances to 7 globular clusters: M5, NGC 6752, M13, M15, 
M92, M30, and M68 obtaining distance moduli of 14.45, 13.17, 14.48, 15.38, 
14.93, 14.95, and 15.29, respectively, with formal uncertainties of 0.1 mag. 
For the last four globulars which are metal poor, these new distance moduli 
are about 0.3 mags higher than the previous standard values. This leads to 
intrinsically bright HB’s (e.g., My = 0.44±0.08 at [Fe/H] = —1.6) with an 
extremely steep My — [Fe/H] slope of 0.5. 

Gratton et ah (1997) have a bigger local sample, but still believe that 
they do not have enough calibrating subdwarfs to make completely model 
independent fits to globular clusters. Instead, using calibrating subdwarfs, 
they determine the relation between [B — V) color for unevolved MS stars 
as a function of metallicity. Then they use this relation to find the loca¬ 
tion of the MS on My - color diagram. Gratton et ah (1997) determine 
distances to 9 globular cluster; 6 coinciding with Reid (1997) [except M15], 
plus 47 Tuc, NGC 288, and NGC 362. Gratton et al. (1997) distances are 
generally slightly shorter than the ones by Reid (1997), but they seem to 
be systematically shorter by 0.2 mag for the three of four metal poor clus¬ 
ters: M30, M68, and M92. Gratton et al. (1997) use the metallicity scale 
of Carreta & Gratton (1997) which is more metal-rich than Zinn & West 
(1984) scale adopted by Reid. All these modifications lead to, 

MviRR) = (0.18 ± 0.09)([Fe/H] ± 1.6) ± (0.45 ± 0.04). (19) 

where we report not the result from original paper but rather the corrected 
relation as reported by Gratton (1998). Equation (19) implies ulmc ~ 
18.60 ±0.07. 
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Studies by Reid (1997) and Gratton et al. (1997) devote a lot of effort 
to correcting for possible systematic effects like Lutz-Kelker corrections, bi¬ 
nary contamination, and errors in reddening determinations. Both studies 
use very similar input data. Nevertheless, they attribute their corrections 
to the pre-Hipparcos MS fitting results to two different sources. Gratton et 
al. (1997) argue that smaller parallaxes to calibrating subdwarfs by them¬ 
selves account for 0.2 mag larger distances to globular clusters. Reid (1997) 
sees the same trend (Hipparcos magnitudes smaller than the ground-based 
measurements), but attributes this change mostly to the lack of sensitivity 
of the sub dwarf luminosity to metallicity for stars with —1.5^ [Fs/H] ~ —2.0 

(which is very unlikely from a theoretical point of view). 

There is, however, some evidence that those two comprehensive studies 
are not the ultimate in MS fitting analysis following the Hipparcos break¬ 
through. Pont et al. (1998) opt for using evolved stars and binaries in their 
analysis, which dramatically increases the number of calibrating stars. They 
find ^m 92 = 14.67, resulting in /? ~ 0.6 [compare to eq. (19)]. Using only 
unevolved calibrating stars, Pinnsoneault (1998, priv. comm.) argues for a 
similar /3 from multicolor analysis of M5 (see also above). 

6.2. WHITE DWARF SEQUENCE FITTING 

The main theoretical difficulty one encounters using WD cooling sequences 
is the dependence of their magnitude on WD mass. Renzini et al. (1996) 
list four independent observations that constrain WD masses: luminosities 
of the red giant branch tip, the HB, the asymptotic giant branch (AGB) 
termination and the post-AGB stars. All of them seem to indicate that 
masses of WD’s in globular clusters (Mcwd) are independent of the host 
cluster metallicity (excellent!) and approximately equal to Mcwd = 0.53 ± 
0.02Mq (Renzini & Fusi Pecci 1988). (Such a narrow range of masses has 
never been either confirmed or falsified observationally, e.g.. Cool, Piotto & 
King 1996 claim that their data for NGC 6397 are consistent with the wider 
range Mcwd = O.55TO.O5M0, although they allow that all the spread may 
originate from the observational errors.) Wood’s (1995) theoretical models 
suggest that the My of the WDs scales with the WD mass as 

5Mv ~ 2A5Mwd (20) 

Therefore, a serious misestimate of the white dwarf My’s requires a very 
severe violation of the WD mass constraints coming from a wide variety 
of theoretical predictions. According to Renzini & Fusi Pecci (1988) the 
My uncertainty produced by differences of WD mass in different globular 
clusters can be of order of 0.05; if a rather conservative range Mcwd = 
0.55 ± O.O5M0 for the mass of the cluster WDs is adopted, 5My can reach 
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0.12. Nevertheless, all observers agree that within a given system, the dis¬ 
tribution of masses is very narrow. This tendency is confirmed by the local 
white dwarfs used to calibrate the My of the WD cooling sequence (Berg¬ 
eron, Suffer, & Liebert 1992; Bragaglia, Renzini, & Bergeron 1995). For 
the solar neighborhood the mean and dispersion of the WD masses are 
(Mwd) = 0.59 ±0.1, but this does not contradict the Renzini & Fusi Pecci 
(1988) claim of a very narrow mass range in globular clusters, where WD 
populations come from low-mass progenitors only. 

There may be also some photometric errors leading to misestimates 
of {B — V), but they are not likely to shift My by more than 0.1 mag 
(in a random direction). Moreover, WDs have practically metal-free atmo¬ 
spheres that come in two varieties: pure H (DA) constituting the great 
majority and pure He (DB), which makes modeling very straight forward. 
They are more abundant than subdwarfs and consequently it is much easier 
to collect a large calibrating sample with accurate trigonometric parallaxes. 
There is, however, one main observational disadvantage of using WD cool¬ 
ing sequences, namely they are at least 10 mag fainter than MS stars of the 
same temperature. For the closest globular clusters, WDs have 23 (De 
Marchi, Paresce, & Romaniello 1995; Richer et al. 1995, 1997; Cool et al. 
1996) and so HST is the best tool to acquire photometric data. 

Even though WD sequences have been observed for at least 6 globular 
clusters (47 Tuc, ui Cen, NGC 6397, M15, M4, NGC 6752), the distance 
to only one globular cluster was determined by the WD fitting technique 
(Renzini et al. 1996 for NGC 6752), one result was mentioned parenthet¬ 
ically (Richer et al. 1997 for M4) and there exist persistent rumors about 
one not yet published (Renzini 1998 for 47 Tuc). Renzini et al. (1996) match 
the WD sequences (when DA is matched, DB falls in the right place) of 
NGC 6752 with the local WD sequence and find hn 6752 = 13.05. When 
combined with the reddening of the cluster and apparent magnitude of the 
HB from Gratton et al. (1997) as originally observed by Penny & Dick¬ 
ens (1986), this yields My{RR) = 0.59 at [Fe/H] = —1.4. Renzini (1998) 
has also recently reported a distance modulus to 47 Tuc that is about 
0.5 mag shorter than the distance obtained through MS fitting and would 
imply My{RR) = 0.95 at [Fe/H] = —0.7. Richer et al. (1997) mention 
that if they use Renzini & Fusi Pecci (1988) cluster white dwarf mass es¬ 
timate of My/D = 0.53 ± O.O2M0, their (dereddened) distance modulus 
to M4 will be 11.09. Adjusting the HB My of Rees (1996) to match the 
Richer et al. (1997) distance modulus to M4, one obtains My = 0.75 at 
[Fe/H] = —1.3. Due to the very small number of analyzed clusters, the WD 
cooling sequence technique must be interpreted cautiosly, but nevertheless 
seems to give some support to the short RR Lyrae distance scale. For the 
3 clusters (and fixing the slope of My - [Fe/H] relation at 0.2), we find 
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Mv{RR) = 0.67 ±0.13 at [Fe/H] = -1.6. 
7. Summary 


TABLE 1. Comparison between different methods to determine RR Lyrae absolute magnitudes 


Method 

Absolute 
magnitude 
at [Fe/H] = -1.6 

Grade 

Main problems 

Potential 

future 

usefulness 

Statistical parallax 

0.77 ±0.13 

A 

— 

A- 

Trigonometric parallax 

0.71 ±0.15 

A- 

done for non-RR Lyrae stars 
small number statistics 

A± 

Cluster kinematics 

0.67 ±0.10 

B± 

modeling of rotation 
uncertain density profile 
proper motions 

A- 

Baade-Wesselink 

0.45-0.70 

B 

temperature scale 
p factor, bolometric corrections 

B 

Theoretical models 

0.45-0.65 

B- 

input physics 

A 

Main sequence htting 

0.45 ± 0.04 

C 

metallicity scale 
reddening uncertainties 
small number of calibrators 

B 

White dwarf fitting 

0.67 ±0.13 

C 

WD masses from theory 
reddening uncertainties 

A 


The most crucial outcome of the Hipparcos mission is the realization 
that astrometry is at the moment the only fully reliable way to measure dis¬ 
tances in the local Universe. However, instead of removing a long-standing 
discrepancy among existing distance scales, Hipparcos results led to numer¬ 
ous analyses which have strengthened old conflicts and created new ones. 
Table 1 summarizes the results reported in this review. Note that the errors 
quoted are statistical only. We have ranked and graded (A - very good; B - 
good; C - acceptable) the methods by our judgement of their susceptibility 
to systematic effects and also included our grade of their potential relia¬ 
bility with forseeable improvement in the data. Of the methods listed only 
3 (statistical parallax, trigonometric parallax and MS fitting) have been 
affected by Hipparcos. 
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